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ABSTRACT 

Aims. We present new data for four candidate obscured Compton-Thick (CT) quasars at z ~ 1-2.5 observed with SINFONI VLT 
spectrograph in AO mode. These sources were selected from a 24/jm Spitzer MIPS survey of the COSMOS field, on the basis of red 
mid-infrared-to-optical and optical-to-near-infrared colours, with the intention of identifying active galactic nuclei (AGNs) in dust 
enshrouded environments, where most of the black hole mass is assembled in dust enshrouded environments. 

Methods. Near infrared spectra were analyzed in order to check for emission line features and to search for broad components in 
the [OIII]-IIj8 and Hq--[NII] regions. X-ray spectral analysis, radio and MIR diagnostics, and SED fitting have also been employed to 
study the nature of the sources. 

Results. We successfully identified three objects for which we had only a photometric redshift estimate. Based on their emission line 
diagnostics and on ancillary multi-wavelength constraints, we find that all four targets harbor obscured AGNs. Broad profiles that 
could be attributed to the effects of outflows are revealed in only one target, MIRO20581. In particular, we clearly resolved a fast 
(~1600 km/s) and extended (~5 kpc) outflow in the [OIII]5007 emission line. This feature, the commonly used indicator for ionised 
outflowing gas, was sampled and detected only for this target; hence, we can not exclude the presence of outflows in the other sources. 
Overall, the constraints we obtain from our targets and from other comparative samples from the literature suggest that these optically 
faint luminous infrared galaxies, hosting obscured AGNs, may represent a brief evolutionary phase between the post-merger starburst 
and the unobscured QSO phases. 

Key words, interstellar medium: jets and outflows - galaxy: evolution - quasars: emission lines 


1. Introduction 


Outflow winds are predicted to be ubiquitous in active galactic 
nuclei (AGN) systems and are invoked in many co-evolutionary 
models to link the growth of supermassive black hole (hereafter 
BH) and galaxies through feedback phenomena. These models 
predict an obscured phase for young recently ignited quasars, 
triggered by the funneling of a large amount of gas into the 
nuclear region during major galaxy mergers (e.g., Menci et al. 
|2008[ [Hopkins et al.|[2008| l. Roughly at the same time, this 
amount of gas is also responsible of vigorous star-formation ac¬ 
tivity. This initial phase is followed by a transitional phase, the 
so-called feedback or blow-out phase ( [Hopkins et al.|[2008] l, in 
which the gas is cleared out through outflowing winds released 
by the BH, before becoming a normal unobscured quasar (QSO). 


* Based on observations with SINFONI VLT spectrograph, ESO pro¬ 
gram 092.A-0884(A) 

** E-mail: michele.pema4@unibo.it 


In the framework described above, during the obscured 
phase, the BH is expected to accrete mass very rapidly, implying 
a vigorous, although obscured. X-ray emission. In the census of 
AGNs, X-ray surveys have been extensively used to probe the 
assembly and growth of BH at high redshift. In particular, since 
the X-ray flux is less attenuated than the optical flux, selection 
criteria based on high X-ray-to-optical flux ratio (/^/ fp) have 
been used to select obscured sources at z ~ 1-2. Several stud¬ 
ies (e.g., Fiore et al.|2003[ Mignoli et al.|2004 Alexander et al.[ 


[2002[[Del Moro et al.|20()9t [Della Ceca et al.|2bl5[ ) have found 

that sources with high fxlfo characterized, on average, by red 
optical-to-near-infrared colours (R-Kvega > 5) and column den¬ 
sities in the X-ray of the order of 10^^“^^cm“^. Moreover, VLT 


XSHO OTER ([Brusa et al.[|20T5l [Perna et al.[[20T5l ) and SIN¬ 
FONI (Cresci et al. 2015j l observations of a small sub-sample 
of obscured QSOs at z~L5, selected on the basis of their ob¬ 
served red colours and high fx! fo ratio, have confirmed the pres¬ 
ence of ionised outflowing material in 75% of object and a dust- 
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R-K (Vega) 


Fig. 1; f 2 A^lml fR as a function of R-K colour for all the COSMOS 
24 pm sources associated with the optical and K-band counter¬ 
parts (green isodensity contours). Black crosses represent the 1 
mJy flux-limited sample. Yellow stars denote the 5 objects ob¬ 


served with SINFONI; cyan stars and squares denote Melbourne 
et al.|(pOTT|l and Brand et al. ( 2007|l targets, respectively (two 


targets present in both samples are indicated with stars); blue 
symbols denote the Perna et ar| ( 2015| l targets. Star and square 
symbols indicate IFU and long-slit observations, respectively. 
Sources marked with red circles are also detected in the 1 mJy 
flux limited sample and are those considered for the analysis. 
The box delimited by vertical solid line and horizontal dashed 
line marks the region proposed in |Fiore et al.| ( |2009) l to select CT 
AGN (Cell A). The horizontal dashed line sets the criterion for 
the selection of DOGS. 


reddened type 1 nature. These objects appear similar to infrared 
bright (Kvega <16) dust-reddened Type 1 QSOs selected by com¬ 
bining radio with near-infrared (NIR) and optical catalogues at 
lower redshift (z < l; |Glikman et al.|[2004| [2007] l. Studying a 
sub-sample of 13 objects, with strongly disturbed morphology, 
[Urrutia et al.| \2Q\2) found that ~60% of these radio-detected 
sources show evidence for outflow in the [OIII]5007 line profile 
(see also|Brusa et al.|2015|l. 


However, each selection technique biases the samples to¬ 
wards particular properties and, chiefly, may biases the charac¬ 
terization of the outflows (see |Brusa et al.|[2015| for details); to 
test co-evolutionary models, we need to select and isolate differ¬ 
ent populations of quasars in the different phases of the AGN- 
galaxy co-evolution, including the initial Compton-Thick (CTQ 
phase. 

Since most of the absorbed AGN energy is re-emitted in 
the mid-infrared (MIR), surveys at these wavelengths can po¬ 
tentially recover the elusive obscured accretion missed by X-ray 
surveys (e.g., Brandt & Alexander|2015|l. Several criteria based 


* Throughout the paper, we distinguish between moderately obscured 
(with log(NH)=22-23 cm“-), highly obscured (with log(NH)=23-24 
cm-2), and CT AGNs (with log(NH) >24 cm^^). 
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on the MIR emission of high-z sources have been introduced 
in the recent years to search for heavily obscured AGNs at z 
~ 1 - 3, and have been applied on Spitzer MIPS observations 
in multiwave length survey fields. Typically, the criteria involve 
the selection of objects with MIR luminosities typical of AGN 


but with faint optical 

or near-infrared emission (e.g., Martinez- 

Sansigre et al.||2005l 

Fiore et al.||20081 |2009t |Dey et al.||2008 

Riguccini et al.|2015|l. 

For example, Fiore et al. ( 

2009]l used the MIPS 2Apm COS- 

MOS catalogue ([Sanders et a 

.|2007|l to select a sample of ~ 60 


candidate obscured AGN/CT QSOs characterized by extreme 
mid-infrared-to-optical flux ratio {fiAum!fn > 1000) in the area 
covered by the C-COSMOS Chandra survey ( |Elvis et al.|2009[ 
|Civano et al.|2012[ ). They coupled this selection with a red colour 
{R - K)vega > 4.5 cut, which corresponds efficiently pick up ob¬ 
jects at the redshift of interest (z ~ 1 - 3). To test the efficiency 
of the selection, they stacked the Chandra images at the position 
of the MIPS sources without a direct X-ray detection and re¬ 
covered a hardness ratio (HR) in the stacking signal larger than 
that measured for less extreme sources if 2 Aiim/fR < 1000 and/or 
(R-K) < 4.5). Still, these sources with f 2 Aiiml Ir > 1000 exhibit 
evidence of both star formation and AGN activity, and there are 
contradicting conclusions on how many of the sources selected 
in this way are actually obscured AGN at z ~ 1-3, rather than 
dusty star forming objects (e.g., [Donley et al.|2008[|Fiore et al.| 
[20091 |Dey et al.|2008^ . 


In this paper we present new SINFONI observations, as¬ 
sisted with Adaptive Optics (AO), for a sample of four luminous, 
highly obscured QSOs in the COSMOS field, selected on the ba¬ 
sis of their high mid-infrared-to-optical flux ratios (MIR/O) and 
red R-K colours. Given the tight positive correlation between 
the f 2 AfimlfR the fxlfo ratios found for sources with col¬ 
umn density of the order of cm“^ (see Fiore et al.|2008 


Fig.2), Fiore et al. (20081 suggested that luminous highly ob¬ 
scured AGNs, that are faint in the X-ray because of the high 
column densities, i.e. N// ^ 10^^“^"^ cm“^, and that cannot be se¬ 
lected using their X/O ratio, can be recovered using their MIR/O 
ratio. Therefore, potentially, in the framework previously de¬ 
scribed, we may be able to select sources in the prelude, or at 
the beginning, of the blow-out phase. 


In the following we will refer to these sources as MIRO tar¬ 
gets, as per their high MIR/O flux ratio, and using their MIPS 


catalogue ID from the COSMOS Spitzer catalogue (Le Floc’h 
jet al.||2009| ). The main aim of the SINFONI observations, be¬ 
sides the spectroscopic determination of the redshifts of the tar¬ 
gets, is to compare the physical properties of the sources selected 
in different ways (MIR/O vs. X/O) and assess the presence of 
ionised outflows and broad features. Targets selected with a sim¬ 
ple MIR/O ratio cut are usually known in the literature as dust 
obscured galaxies (DOGs, |Dey et al.|2008| see also Riguccinretj 
|aL]|201 l||2015| for a complete discussion on the DOGs popula¬ 
tion in COSMOS). We will compare our results with a compila¬ 
tion of few DOGs sharing similar properties, e.g., redshift, 2Apm 
flux, in the last Section. 


The paper is organised as follows: Sect. 2 presents the sam¬ 
ple selection and the ancillary data collected for our MIRO tar¬ 
gets; Sect. 3 outlines the VLT observations and data reduction; 
Sect. 4 exposes the spectroscopic analysis. Sect. 5 presents proof 
of ionised outflowing material in the X-ray source MIRO20581 
and discusses the energetic output associated to the outflow and 
Anally we summarise our results and the implications in Sect. 
6. Throughout the paper, we adopt the cosmological parame¬ 
ters //o=70km/s, Qm= 0.3 and 12^=0.7 ( jSperger et al.|2003] l. We 
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Table 1: MIRO/SINFONI sample: selection properties; log file of observations 


MIRO 

(1) 

RA 

(2) 

DEC 

(3) 

^phot 

(4) 

R 

AB 

(5) 

K 

AB 

(6) 

f24fim 

mJy 

(7) 

R-K 

Vega 

(8) 

MIPS/O 

(9) 

band 

(10) 

Guide Star 
(name) 

(11) 

expo 

(min) 

(12) 

Zspec 

(this work) 
(13) 

18744 

10:01:52.2 

01:56:08.6 

0.97“ 

22.8 

19.6 

2.08±0.03 

4.90 

780 

J 

Hip037044 

50 

0.97 

10561 

09:59:43.5 

01:44:07.6 

1.54 

24.6 

21.4 

1.70±0.06 

4.80 

3090 

J,H 

Hip040661 

50,50 

1.43 

28704 

10:01:45.9 

02:28:53.8 

1.74 

26.3 

21.6 

1.63±0.02 

6.35 

14800 

HK 

Hip046054 

30 

1.64 

20581 

10:00:00.6 

02:15:31.1 

2.09 

25.3 

22.6 

1.47±0.02 

4.35 

5180 

HK 

Hip046054 

80 

2.45 

18433 

10:01:44.8 

01:55:55.8 

2.59 

24.6 

21.6 

0.44±0.02 

4.75 

880 

HK 

Hip044598 

120 

— 


Notes: (1) target name; (2) right ascension; (3) declination; (4) photometric redshifts available prior to the SINFONI observations; (5) and (6): R 
and K-band magnitudes; (7) MIPS lAjim flux; (8) R-K colour; (9) MIPS 24pm/0 flux ratio; (10) SINFONI filters; (11) guide star name; (12) total 
integration time on-target for each band; (13) spectroscopic redshift. 

“ Spectroscopic redshift available from Magellan IMACS spectrum. 


Table 2: MIRO/SINFONI sample: main properties 


MIRO 

(1) 

Nh 

102'’ cm-2 
(2) 

log(L2_io) 
erg s“' 

(3) 

l0g(L5 8) 
erg 
(4) 

log(Lb„i) 
erg s“* 

(5) 

Sradio 

ply 

(6) 

924 

(7) 

M. 

1O'‘M0 

(8) 

SFR 

Mo/yr 

(9) 

log(R,B)'’ 

(10) 

E(B-V)a,„ 

(11) 

E(B-V)^Civ 

(12) 

18744 

2.4!'-= 

43.82 

44.85 

45.11 

424±28 

2.3 

2.0 

99 

0.2 

0.5 

2.7 

10561 

2.9+JI 

43.96 

45.35 

45.97 

72±15 

0.4 

1.5 

196 

0.3 

0.7 

4.0 

28704 

— 

— 

45.32 

47.63 

154±25 

2.0 

0.6 

11(<25)‘' 

-0.8(<-0.5)'’ 

0.5 

4.1 

20581 


45.00 

45.89 

46.61 

5430 ±60 

-0.6 

1.9 

48(<132)“ 

-0.6(<-0.2)“ 

0.6 

2.0 


Notes: (1) target name; (2) column density; (3) intrinsic X-ray luminosity; (4) rest-frame 5.8yum luminosity; (5) bolometric luminosity; (6) rest- 


frame 1.4 GHz flux; (7) q 24 oij = log(f24p,ml fiACHz) (Bonzini et al.|2013|l; (8) stellar mass; (9) star formation rate; (10) starburstiness, defined as the 
ratio between the specific Star Formation Rate (sSFR=SFR/M.) and that expected for Main Sequence galaxies at given computed stellar mass and 


spectroscopic redshift (sSFRms,z), according to the relation ofjWhitaker et al.H2012f; (11) galaxy reddening ; (12) AGN reddening. 
“ Values in the parenthesis refer to measurements constrained using upper liniits in the FIR SED (see Figure]^. 


adopt a Chabrier initial mass function to derive stellar masses 
and star formation rates (SFRs). 


2. Sample selection 

Figure [T] shows the mid-infrared-to-optical flux ratio f 2 AumlI r 
versus me (R-K) colour diagnostics diagram proposed in |Fiore et 
al. ( 2008|l, applied to th e MIPS-selected sources in the COSMOS 
field ( Fiore et al.|2009[ ). The green isodensity contours show the 
distribution of the full sample of MIPS-selected sources in the 
COSMOS G03 data ( |Le Floc’h et al.||2009| l and with associ¬ 
ated optical and K-band counterparts (~ 15000 sources), while 
the black crosses are the sources detected at fluxes larger than 
1 mJy at 24pm (~ 550). The box in the top right corner in the 
colour-colour space, delimited by vertical solid line and horizon¬ 
tal dashed line, mark the region where CT AGN are expected in 
more than 60% of the MIPS-selected sources (cell A in Fiore et 
al. 2009). All the sources above the horizontal dashed line are 
instead usually referred as DOGs in the literature. 

From the |Le Floc’h et alT ( 2009) 1 MIPS selected sample, we 
observed with SINFONI-NGS mode five targets marked as yel¬ 
low stars in Figure [1] and listed in Table [1] with the fluxes in 
the R, K and MIPS 24pm band used for the selection. Three 
out of five are classified as DOGs (MIRO10561, MIRO28704, 
MIRO20581) and 2 out of five are also in Cell A of Fiore et al. 
(2009; MIRO10561, MIRO28704f] The main selection criteria, 
in addition to the high MIPS/O fluxes and the red R-K colour, 
were the proximity to a bright AO star and a photometric red¬ 
shift broadly in the range z~ 1 - 3 so that rest-frame optical lines 


^ At the time of observations, all the targets were selected within the 
cell A. Differences from the current situation shown in Figure are 
due to the fact that we now use an improved version of the COSMOS 
photometric catalogue (Laigle et al. 2015 in preparation) in which the 
photometry is slightly changed. The colours in the figure are all related 
to the “total flux” measurements. 


ACS-HST I 


r+ 


Ks 



IRAC3.6um MIPS 24um 


@ 

® 


5 ^ 


Fig. 2: From left to right; ACS-HST I, r+ SUBARU, Ks COS¬ 
MOS, IRAC 3.6 pm, MIPS 24 pm band cutouts of the five MIRO 
targets. The target name, the position, the cutouts scale and the 
redshift derived from the SINFONI data are also labelled. For 
display purposes, 3"x3"ACS-HST cutouts show the regions in 
the red boxes superimposed on the r-t cutouts. 


are redshifted in the SINFONI J, H or K filters (four targets). 
In addition, for only one source (MIRO 18744) a spectroscopic 
redshift was available from the IMACS/Magellan follow-up of 
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X-ray sources in the COSMOS field (Trump et al. 2007), and 
we proposed to observe its Ha'+[NII] region with SINFONI (J 
band). 

In Figure [T we also show the z ~ 2 DOGs presented in 
Melbourne et al. (2011; cyan stars) and in Brand et al. (2007; 
cyan squares) for which K-band magnitude is available, and the 
z ~ 1.5 dust-reddened type 1 sources presented in Perna et al. 
(2015; blue symbols). The properties of these targets will be dis¬ 
cussed in Section in order to compare our results with those 
previously reported in literature. 

The ACS (3"x 3"), R-band, K-band, IRAC3.6pm and MIPS 
(15"xl5") cutout^of the 5 targets are shown in Figure ^ We 
included also the IRAC 3.6/tm cutout in order to verify blend¬ 
ing problems in the MIPS emission; for all but MIR018433, 
we can safely say that the majority of the emission at 24 mi¬ 
cron is correctly associated to the K-band and optical counterpart 
(e.g. our SINFONI targets), and that the observed colours do not 
suffer from bad photometry. MIRO18744 may show evidence 
for an ongoing merger: tidal tails and double nuclei are distin¬ 
guishable in the ACS cutout. MIR018433, instead, presents two 
components in the ACS cutout (last row of Figure |^, which 
are strongly blended at optical and infrared wavelengths and 
preclude an accurate SED fitting decomposition and a correct 
photometric redshift derivation. Indeed, MIR018433 is the only 
source for which no spectral features have been detected in the 
SINFONI HKband in 2 hours observations (see Table[T]i. More¬ 
over, MIR018433 is the only source below the ImJy flux lim¬ 
ited sample (red circles in Figure 1). For all the above reasons, 
MIR018433 was excluded from the subsequent analysis. 


2.1. Identikit via ancillary data 

Previous works (e.g. Alexander et al. 2002, Donley et al. 
2008) have demonstrated that red optical-to-near-infrared 
colours and high MIR-optical ratios can identify both AGN and 
star-forming galaxies. In the following, we briefly discuss the 
multiwavelength properties of the SINFONI targets in order to 
assess which one among SF or AGN activity is the dominating 
process. We note however that all our four targets have / 24 ;j,„ > 
1 mJy, and that the 24/im emission is, on average, increasingly 
dominated by AGN contribution at higher f 24 /jm (s-g-. Brand 


|et al.|[200^ |Dey et ar]|2008| l. In the subsequent analysis (e.g. 
X-ray and SED fits) we made use of the spectroscopic redshift 
obtained from our SINEONI observations. 


1. MIR-X-ray diagnostics: The AGN intrinsic hard X-ray lumi¬ 
nosity and the infrared luminosity re-emitted by the torus fol¬ 
low a tight correlation ( |Lutz et al.|2004[ [Gandhi et al.|2009| l. 
Eigure shows the distribution of rest-frame observed X- 
ray luminosity (L^*') vs. for several sampl es of CT 

candidate AGbf] collected by|Lanzuisi et al.||2015a (with the 
addition of XMM ID 5371; |avano et al.|2()15| ) md for our 
MIRO targets (red diamonds). 

Eor all the sources, the rest-frame 5.Spin luminosities were 
obtained using a simple power-law interpolation between the 
24 and Spm observed-frame luminosities. The X-ray lumi¬ 
nosities have been computed on the basis of the available 


^ The cutouts are extracted from the public COSMOS cutouts 
page: http://irsa.ipac.caltech.edu/data/COSMOS/index_cutouts.html 
All the sources having upper limits to the X-ray luminosity have been 
pre-selected as AGN candidate using a variety of methods specifically 
designed to discriminate between SF and AGN galaxies. 
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Eig. 3: Log(Lj*^ ) vs Log/Ls.g^m) for several CT candi¬ 
dates. Grey circles represent literature CT candidate AGN (see 
[Lanzuisi et al.|p015a| for more details). Red diamonds repre¬ 
sent our MIRO targets; intrinsic X-ray luminosities of the X-ray 
detected objects are also indicated with upper bars. The green 
shaded area is the relation of Lutz et al. (2004) for a sample of 
low redshift unobscured AGN. The orange solid line is the rela¬ 
tion for high redshift unobscured AGN (Eiore et al. 2009), while 
the dashed and long dashed lines are the expected relation for a 
10 ^^ cm“^ and 10 ^^ cm“^ absorber. 


XMM and Chandra data in the COSMOS field. More in de¬ 
tails, MIRO20581 and MIR018744 are both detected in the 
X-rays, in the XMM-COSMOS (XID70135 and XID60205: 
ICappelluti et'ai:||2007| |Brusa et al.||20T0l l and C-COSMOS 
(CID451 and CID401; Civano et al. 2012) surveys, while 
MIRO10561 is detected in the deeper COSMOS-Legacy 
survey (CID 3587 in the Civano et al. 2015 catalogue 
and Marches! et al. in preparation). The remaining source, 
MIRO28704, is instead undetected down to a luminosity of 
log(Lx)~ 42.8 in the 2-10 keV band. 

The Lutz et al. (2004; green shaded area) and Eiore et al. 
(2009; orange solid line) relations represent the tight X-ray- 
to-mid-IR correlations found for the low- and high-redshift 
(z>l) unobscured AGN, respectively. The Lutz et al. (2004) 
relation has been confirmed recently also at higher redshift 
by |Mateos et al.|p015| l. These relations, however, have been 
calibrated at low luminosities, and require extrapolations to 
high luminosities (i.e. logiAL^iS.Spin)) > 46). A flattening 
of the MIR-Xray relation at the highest luminosities has been 
found by the recent work of Stem| ( |2015[ ). Assuming that 
both the hard X-ray and the infrared luminosities are re¬ 
lated to the AGN activity (see below; Eigure [^, given that 
the mid-IR is largely independent of obscuration, a lower 
Lx to L 5 8 ratio with respect to that observed for unobscured 
AGN suggests that the observed Lx are affected by obscura¬ 
tion (e.g. Stern et al.|2014 l. All the MIRO sources lie below 
these relations, at values consistent with a heavily obscured 
absorber; the long- (short-) dashed line in Eigure[^marks the 
Nh = 10^“* (10^^) cm“^ locus, computed from the Eiore et al. 
relation. Therefore, the Nh loci in the figure are computed in 
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Energy (keV) 


Energy (keV) 
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Fig. 4: X-ray spectra of MIRO20581 (left), MIR018744 (centre) from XMM and Chandra, and MIRO10561 (right) from Chandra. 
The XMM (black) and Chandra (red) data of MIRO20581 and MIRO18744 are fitted with a double power-law (as shown in the 
model components). The MIRO10561 Chandra data are instead fit with a single absorbed power-law. 


the most conservative approach: using the Lutz et al. relation, 
the dashed lines would be steeper than those obtained from 
the Fiore et al. one. Therefore, the intrinsic X-ray luminos¬ 
ity would be larger and, as a consequence, the Nh needed to 
explain the observed luminosities of the MIRO targets would 
also be larger. 

2. X-ray spectra: The X-ray spectra of the three X-ray de¬ 
tected MIR/O are shown in Fig. Given the low photon 
statistics available for all the detected sources (in the range 
20-150 net counts) we applied the Cstat statistic ( |Cash|1979] l 
to the unbinned data and assumed a very simple model in 
order to recover a rough estimate of the nuclear obscura¬ 
tion and instrinsic luminosity: a power-law with photon in¬ 
dex fixed to F = 1.9 plus obscuration at the source redshift 
(plus galactic Nh). 

In the XMM-detected sources, a second component is 
required by the data in order to model the soft emission. 
Indeed, as demonstrated in |Lanzuisi et al.| ( |2015a| l, given 
the complexity of the X-ray emission, and the concurrent 
presence of other processes such as scattering compo¬ 
nents or emission from star formation, heavily obscured 
AGN can be missed when fitting low-counting statistics 
data compatible with a single power-law model (see also 
Lanzuisi et al. 2015b). For all these sources we derived 
column density of Nh ~ 2.5-7x 10^^ cm“^, with high 
uncertainties (see Table 2). The X-ray detected sources, 
therefore, although characterized by a high obscuration, are 
not in the Compton Thick regime, as expected given current 
X-ray surveys limits and sensitivities and consistent with 
previous works (e.g., [Lanzuisi et al.|2009[|Georgakakis et al.| 
2010| l. The rest-frame intrinsic X-ray luminosities are also 
listed in Table |2l The most luminous of the three sources 
is MIRO20581 with an X-ray luminosity ~ 10"^^ erg s“\ 
while the other 2 sources have inferred intrinsic luminosities 
slightly below 10"^^ erg s“'. These luminosities exceed of 
about 2 dex those expected from stellar processes given the 
observed SFR (see below). 


3. SED fitting: Figure shows the SED fitting decomposition 
of the four targets considered in the paper, obtained mak¬ 
ing use of a modified version of the magphys code ( |da| 
|Cunha et al.||2008] ) designed to take into account a possible 
AGN emission component ( Berta et al.|2013 1 together with 
a modestly-absorbed galaxy component. They have stellar 
mass M, in the range 0.6-2x10" Mq and SFR of 10-100 
Mo/yr, and are in (or below) the main sequence (MS) of star 
forming galaxies (see [Whitaker et al.|2012|l, as suggested by 


their ’starbustiness’ R 5 £=sSFR/sSRFa/s_ 5 , (see Table |^. For 
two sources, MIRO20581 and MIRO28704, the far-infrared 
emission is not well constrained (see Figure]^, hence in the 
table are also reported SFR and R^ b upper limits in parenthe¬ 
sis, computed treating the FIR upper limits as real detections. 


In all cases the (observer frame) NIR emission is dominated 
by the host galaxy light and the AGN disk component is con¬ 
siderably extincted (E(B-V)=2-4). According to the SED fit¬ 
ting decomposition, the 5.8/im luminosity is dominated by 
the torus emission = 87%), in agreement 


with the results from the bright DOG sample presented in 
[Riguccini et al.| ( 2015| l. This also confirm that the X-ray-to- 
mid-IR diagnostic discussed above is a reliable instrument to 
test the X-ray obscuration of the MIRO targets. 


According to the criteria exposed in Dey et al. (2008i 
who classified dust obscured galaxies in AGN-dominated 
(“power-law” DOGs) and SE-dominated (“bump DOGs”) on 
the basis of the rest-frame optical-to-mid-infrared SED shape 
(see their Section 3.1.2, and their Eig. 5), all the four MIRO 
targets would be classified as power-law DOGs, although 
MIRO 18744 appear to show intermediate characteristic be¬ 
tween the two classes. 


4. Radio: All the targets are also detected in the Very Large Ar¬ 
ray (VLA) observations of the COSMOS field ( [S chinned 
[et al.|2010) l. [Fiore et al.[ ( |T009[ ) reported that QSOs selected 
on the basis of the MIR/O excess at z~L5 are more ra¬ 
dio luminous than unobscured type 1 QSOs of similar lu¬ 
minosity and redshift, when the intrinsic 5.8 pm luminosi¬ 
ties are compared (see also [Martmez-Sansigre et aT][2005| l. 
MIRO20581 is the only radio loud target {q 2 Aohs — -0.6, 
being q 24 obs = log{f2Aiimlf\.ACHz)\ see [Bonzini et al.[[2013 


Fig.2) in our SINFONI sample, with Li.4Gft=4.5 x 10^ W 
Hz '. MIRO20581 is also detected in the 3GHz survey of the 
COSMOS field (Smolcic et al. 2015, in prep) and is one of 
the most luminous sources in the Very Long Baseline Array 
(VLBA) COSMOS catalogue (Herrera-Ruiz et al. 2015, in 
prep); its Jansky Very Large Array (JVLA) measurements at 
3 and 1.4GHz imply an inverted radio spectral index, which 
is consistent with a compact radio source rather than with 
a diffuse star forming region. Indeed, using the relation in¬ 
troduced by Condon (1992 1 , between the 1.4 GHz luminos¬ 
ity and the SFR, we derive for MIRO20581 a SFR a fac¬ 
tor of 10-40 larger than the SED fitting estimate, and there¬ 
fore the radio luminosity is interpreted as due to AGN ac¬ 
tivity. MIR018744 and MIRO28704 are also detected in the 
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VLBA catalogue and even in these cases, possible signatures 
of compact radio cores are present. 

5. High ionisation diagnostics'. For MIRO18744 a IMACS 
Magellan spectrum is available ( |Trump et al.||2007] ), sam¬ 
pling the rest-frame range 2850-4600 A. The spectrum 
shows a prominent [NeV]3425 emission line, an unambigu¬ 
ous sign of obscured nuclear activity ( [Mignoli et al.|[20T3j 
[Lanzuisi et al.|2015b| ). 

The main properties of our SINFONI sample discussed 
above are reported in Table with the targets sorted by decreas¬ 
ing MIPS flux and increasing redshift. Overall, the multiwave¬ 
length constraints we have in our lApm bright SINFONI targets 
suggest unambiguously the presence of obscured AON activity. 
These characteristics correspond to those expected for objects 
caught in the post-merger dust-enshrouded phase of rapid black 
hole growth (see e.g. |Hopkins et al.||2008l |Fiore et al.||2008| 
[2009] !. 


3. SINFONI observations and data reduction 


The observations were obtained in service mode using the near- 
infrared spectrometer SINFONI of the VLT in adaptive-optics 
(AO)-assisted mode, during period 92A (from 2013-12-28 to 
2014-03-30). All the targets in the sample were observed in one 
or two of the SINFONI Alters (J, H, K or HK) depending on the 
initial redshift guess. We note that the program has been exe¬ 
cuted only partially (~ 45%) and therefore our targets have not 
been observed with all the requested Alters and/or for the entire 
requested time. We used a fleld-of-view (FoV) of 8x8” in a 2D 
64x64 spaxel frame. The spectral resolutions are R~1800 for J, 
R~2900 for H and R~1400 for HK. 

We achieved a spatial resolution of 0.2" (FWHM) based on 
the point spread functions (PSF) obtained in natural guide star 
(NGS) AO-mode, which roughly corresponds to 0.9 kpc at the 
average redshifts of z=1.5. This spatial resolution is in agree¬ 
ment with those obtained in other SINFONI AO-assisted obser¬ 
vations (e.g., |Bouche et al.|2013[|Cresci et al.|2009| l. Our targets 
do not extend more than ~ 1-2” in diameter, and were there¬ 
fore observed with on-source dithering in order to use the object 
exposure with the closest MID as an approximation of a sky ex¬ 
posure. The information about the observations for each object 
are shown in Table [T] 

Besides the objects of the sample, a set of standard stars and 
their respective sky frames were also observed to flux-calibrate 
the data. Guide star names are also reported Table 1. The stars 
have R magnitudes in the range 15<R<16. 

The data reduction process was performed using ESOREX 
(version 2.0.5). We used the IDE routine "skysub.pro” ( [Davies | 
|[2007] l to remove the background sky emission. Then, we used 
our own IDE routines to perform the flux calibration and to re¬ 
construct a final data cube for each object, adding the different 
pointings. The flux calibration was performed following the pre¬ 
scription indicated byjPiqueras L6 pez et al. (2012|l. 


4. Data analysis and Spectral fits 

Here, we briefly discuss the general data analysis and results of 
the spectral fits. In the next sections we describe the more de¬ 
tailed spatially-resolved analysis for our best-case, MIRO20581. 

Eigure shows the one-dimensional integrated spectra, ex¬ 
tracted in a 1-2” diameter aperture, according to the compact¬ 
ness of the source. All the targets show the HaH-[NII] complex, 
with the exception of MIR0287()4 which has a lower-quality 
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spectrum and in which only the Hci[^ emission is detected. The 
only target for which we detected also the [OIII]5007 emission 
line is MIRO20581. MIRO10561 was observed in the J band 
too, but both [GUI] line and continuum emission were not de¬ 
tected. Instead, [OIII] emission line for the other two targets was 
not covered, because of the incompleteness of the observations 
(MIRO28704; see Section and of the wavelength coverage 
of SINEONI instrument (MIR018744; [OIII]5007 expected at 
'^obs-frame ^ 0.986///72). 

We fitted simultaneously each of the emission lines (from 
only one, Ha in MIRO28704 to a total of six, H/3, [OIII] doublet, 
[Nil] doublet and Ha, for MIRO20581) with Gaussian line pro¬ 
files. When more than one emission line is fitted, we constrained 
the centroids and the line flux ratios according to atomic physics, 
while the widths were fixed to be the same as in each emission 
lines (see jPerna et al.|2015| l. 

From the fit described above, we computed spectroscopic 
redshifts for our targets. We chose as redshift solution the 
one which produces the best fit of the narrow components of 
the emission lines. We detected Ha and [Nil] emission for 
MIRO 18744 at z=0.97, consistent with the spectroscopic red¬ 
shift already available from the Magellan spectrum. For the re¬ 
maining three sources, we were able to assign for the first time 
a spectroscopic redshift from our line fit (see column 13 in Ta¬ 
ble!^. The spectroscopic redshifts are in general agreement with 
the photometric estimates available within the COSMOS sur¬ 
vey ( jSalvato et al.||20lT| jllbert et ^|2009| ), with an accuracy 


ttf \Zphot ^5/7ecl/(l Zspec) ^ 0.1. 

The results of the emission line fits are reported in Table 
1^ In order to investigate the nature of the ionising source, we 
investigate the emission line ratios diagnostics [NII]/Ha and 
[OIII]/Hy8. The only source for which we have both [NII]/Ha and 
[OIII]/H/3 ( MIRO20581) to calculate the BPT diagram ( jBaldwinj 
jet al.|198T l, lies in the AGN photoionisation region in such di¬ 
agnostic. For the remaining sources, although with large uncer¬ 
tainties due to the low-quality spectra, two out of three also show 
[NII]/Ha ratios consistent with an AGN origin ( jVeilleux & Os-j 
|terbrock|1987| ), consistent with the AGN classification discussed 
in Section 2.1. The diagnostic line ratios are reported in Table[^ 
lower limit in the [OIII]/HjS ratios are due to the non-detected 
HyS emission line. 

We note here that although the clear presence of two objects 
in the ACS-HST image of MIR018744 (Figure]^, given the low 
S/N, the SINFONI spectrum is relative to both objects. 


5. The ionised outflow in l\/IIRO20581 


We now concentrate on MIRO20581. In order to reproduce the 
line profiles in the [OIII] and Ha region we had to introduce an 
extra blueshifted and broad (FWHM^ 1600 km/s) component 
(see Figure]^ panel a). Emission from ionised gas in forbidden 
lines like [OIII] can not be associated to motion in the broad line 
region (BLR) because it would be suppressed by collisional de¬ 
excitation when produced in high-density regions (but see jPev-j 
ereux|20iT for an alternative explanation). Eor this reason, any 


^ The narrow feature at A «6584 A is associated with a wrong sky¬ 
line subtraction. We discard the possibility that the observed line is 
HeI10830 at z=0.6. We fitted the SED imposing this redshift. This fit 
produced a significantly larger Chi-square, 4 times the value obtained 
imposing z=1.64. Indeed, at z=0.6 the source would be undetected 
down to a luminosity of log(Lx )~ 42.2 in the 2-10 keV band; with 
a bolometric luminosity estimated by this SED fitting decomposition 
of log(Lj„/)=45.3, MIRO28704 would have a very unusual bolometric 
correction of the order of k;,„;=Li,o//Lx > 1000. 
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M. Pema et al.: SINFONI spectra of heavily obscured AGNs in COSMOS: evidence of outflows in a MIR/0 target at z~ 2.5 




Fig. 5: SED fitting of MIRO targets. The black dots represent the observed data points (in each panel, from left to right: SUBARU U, 
B, g, V, r, i, z; VISTA Y, J, H and Ks; the four Spitzer IRAC bands, Spitzer MIPS24; Hershel PACS and SPIRE). The blue line shows 
the integrated extincted emission originating from the host galaxy. The magenta line represents the star formation contribution for 
dust absorption, partially redistributed across the MIR/EIR range in a self-consistent way (da Cunha et al.|2008 Berta et al.|2013 


Delvecchio et al. 2014|l. The red line reproduces the AGN contribution and incorporates both the accretion disc and the torus 


emission. The black solid line represents the sum of all components. 


broad (> 550 km/s) profile in such forbidden lines is generally 
interpreted as to be ascribed to outflowing ionised gas. Broad 
Ha profiles may be also ascribed to the presence of a ionised 
outflowing gas (e.g. |Genzel et al.|2014| l; however, the Ha emis¬ 
sion may suffer from severe contamination by the presence of the 
BLR motion and therefore may be considered less reliable tracer 
of outflows and the associated energetics. This is especially the 
case when the high-velocity BLR wings are not detected, as in¬ 
stead, for instance, it is the case of Mrk 231 (|Rupke & Veilleux| 

[20TTT l. 

In order to investigate at best the ionised gas emission in 
MIRO20581, and taking into account the considerations above, 
we fitted simultaneously the H/l-i-[OIII] and Ha-i-[NII] regions 
(see Section]^ for details) using two approaches. In order to re¬ 
produce the line profiles of all the emission lines, we fitted the 
two regions 

1. with two sets of Gaussian profiles: one to account for the 

presence of NLR components (in the following narrow com¬ 


ponent, NC), with EWHM<550 km/s, one for the presence 
of outflow components (OC), with EWHM>550 km/s; 

2. using the same components as above, namely the NC and OC 
components, and adding a broader profile ( EWHM>1900 
km/s) to account for the presence of the Ha emission origi¬ 
nated in the BLR. 

The two fits have acceptable low residuals and adequately 
represent the shape of the line profiles, hence it is not possible 
to confirm or exclude the presence of the BLR emission in this 
source, at least with this S/N. The results of these two approaches 
are reported in Table If present, the Ha BLR emission would 
have a EWHM of about 2500-3000 km/s. The lower value has 
been obtained fitting the entire Ha-i-[NII] profile with only one 
Gaussian. The NII/Ha ratio depends on the detailed modeling 
but, in both cases, the emission lines remain in the AGN photo¬ 
ionisation region (see previous Section). 

A tentative fit with 2 sets of Gaussian profiles, taking into 
account only NLR and BLR components was performed. This 
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Fig. 6; (a) MIRO20581 HKband spectrum around the [OIII] (left) and the Ha+[NII] complex (right). Superimposed on the spectrum 
are the best fit components (solid and dashed blue curves, with arbitrary normalization in order to ease the visualization). The red 
solid curves represent the sum of all components, including the power-law. Dotted lines mark the wavelengths of the HyS, [Oil] 
doublet, Ha, [Nil] and [SII] doublet, (b) MIRO10561 H band integrated spectrum around the Ha-i-[NII] complex.(c) MIR018744 
J band integrated spectrum around the Ha-i-[NII] complex.; (d) MIRO28704 HK band integrated spectrum around the Ha-i-[NII] 
complex. See (a) for the description of superimposed curves in (b), (c) and (d) panels. In the bottom panel of each fit the residuals 
with respect to the best fit are shown. 


has shown high residuals especially in the [OIII]5007 profile. 
Figure [7] shows these residuals. 

5.1. Spatial analysis 

Figure panel a, shows the contour plot of the median SIN- 
FONI data-cube over the entire HK wavelength range, in steps 
of Icr starting from 3cr (green scale). The standard deviation cr 
was computed in a 1.25-1.75” annulus centred on the target. The 
astrometry in the SINFONI data-cube was performed using the 
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Ks COSMOS cutout, obtaining a match between the coordinates 
of the peak of intensity in the Ks COSMOS cutout and in the me¬ 
dian SINFONI data-cube. In Figure 7a, the Ks COSMOS con¬ 
tours are over-imposed on the SINFONI data-cube (cyan solid 
curves). We also show the HST/ACS F814W contours at higher 
resolution (magenta curves; see also Figure 1). 

Figure panel b, shows the [OIII]5007 map integrated on 
the continuum-subtracted total line profile (grey scale), with 
over-imposed contour levels (starting from 3 cr) of the emission 
of the line core (5000-5014A; width ~ 800 km/s) and of the 
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M. Pema et al.: SINFONI spectra of heavily obscured AGNs in COSMOS: evidence of outflows in a MIR/0 target at z~ 2.5 


Table 3: Emission line properties in the integrated spectra 


MIRO 


ha 

flOHI] 

fua) 

log(/[o///]/Hy8) 

FWHM 



(10 erg/s/cm^) 



(km/s) 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

18744 


32.4±4.9 

- 

0.12±0.08 

- 

565±20 

10561 


11.8±1.3 

- 

-0.10±0.07 

- 

560±25 

28704 


9.4±1.9 

- 

- 

- 

430±100 

20581 

NC 

4+3 

^-2 

9 c+2.5 
^•‘^-1.5 

0.88/“iS 

> -0.47 

360±70 


OC 

43.5±14.8 

95±6 

-0.16±0.51 

> 0.47 

1600±100'' 


NC“ 

4+3 

^-2 

9 c-i-2.5 
^•‘^-1.5 


> -0.47 

370±75 


OC" 

24±13 

95±6 

-0.28±0.6 

> 0.47 

1600±100'' 


Notes: (1) target name; (2) kinemetic component: NC = narrow component, OC = outflow component; (3) Ha flux; (4) [OIIIJ5007 flux; (5) 
diagnostic [NII]6583/Ha flux ratio; (6) diagnostic [OIII]5007/H/1 flux ratio; (7) kinematic component width. 

“ These results were obtained following the second approach illustrated in Sec.|^ adding a BLR component for the Ha profile. 

* The centroid of the OC profile is blueshifted of ~ 550 km/s. 
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Fig. 7; MIRO20581 residuals in the [OIII] and Halpha range 
obtained by fitting the observed lines with only NLR and BLR 
components (See Section 5; to be compared to Panel a of Figure 

E- 


blueward part (4973-5000A; width~ 1600 km/s) of the line pro¬ 
file. Blueshifted emission is found out to a distance of R=0.6” 
(associated to region B), i.e. 4.8 kpc from the nucleus. Instead, 
the core emission is less extended, and could be associated with 
the narrow component of the [OIII] profile. 


5.2. Kinematic anaiysis 


In order to map the line emission distributions and the corre¬ 
sponding velocities, a nuclear and a off-nuclear spectra are ex¬ 
tracted from two 3x2 spaxel regions (3x2 kpc) close to the central 
and the off-nuclear peaks in the [OIII] channel map ( Fig.[8|panel 
b, regions A and B labelled with black boxes; see also |Pema et al.| 
2015[ Cresci et al.|2015| l. In Figure]^ panel (c), we show the in¬ 
tegrated spectra over the two regions: red and orange solid lines 
represent the nuclear and off-nuclear line profile respectively. 

In both cases, the [OIII] profiles are broad, with widths of 
FWHM= 1400-1200 km/s, and maximum velocities of 1600- 
1650 km/s (nuclear-offnuclear spectrum). In the off-nuclear re¬ 
gion, the centroid of the emission lines is blueshifted of ~700 
km/s. 


5.3. Outfiow properties 

Assuming that the broad and shifted [OIII] component can be 
associated with an outflowing wind, the kinetic power (P^") and 
mass-outflow rate of the outflow can be computed under 


reasonable assumptions in the case of a biconical geometry. First 
of all, given that the electron density of the outflowing gas can 
not be estimated directly from the data due to the low quality of 
the spectra in the [SII] region (see Figure]^, and the metallicity 
indicators ( |Pettini & Pagel||^04[ |Yin et al.]|2007] l are not use¬ 
ful because of the AGN ionising radiation, we assumed standard 
values of (100 cm“^) and metallicity (solar). Generally, the 
HyS luminosity is used to calculate the total amount of gas and the 
mass outflow rate (e.g., Liu et al.|2013[[&esci et al.|2015 1. How¬ 
ever, in our case this line is not detected, and we have adopted 
the |Cano-Dfaz et 3T] ( |2012| l formulae which instead employ the 
[OIII] line luminosity: 


CL44([0///])t^L,3 
7 X 10 - ^rg s 

nesRkpclOiom ^ 


( 1 ) 


M'Z = 164 X 10 


43 CL44([OIII])Vout,3 

rie^RkpclO^om 


Mq s' 


-1 


( 2 ) 


where L 44 ([Oin]) is the [OIII] luminosity associated to the 
outflow component in units of 10"^"^ erg s \ is the electron 
density in units of 1000 cm“^, Vout ,3 is the outflow velocity Vgut in 
unit of 1000 km s“\ C is the condensation factor (:« 1), lO^'^^^i 
is the metallicity in solar units, R^pc is the radius of the outflow¬ 
ing region in units of kpc. We therefore used the [OIII]5007 flux 
associated to the outflow component in the 1” integrated spec¬ 
trum (Figure]^ panel a). We further adopted a spatial extension 
of 4.8 kpc for the outflowing gas given that we observe the blue- 
ward emission out to this distance (see Figure]^ panel section 
5.1). Finally, we considered as outflow velocity the maximum 
velocity observed Vma.x in the nuclear region {vout = 1600 km/s; 
see Section |5.21 i, and we assumed that lower velocities are due 
to projection effects ( |Cano-Dfaz et al.|2012[|Cresci et al.|2015| l. 

Following the Cano-Dfaz et al. formalism, the kinetic power 
is P'°''-l.5x 10“*"^ erg s“', while the outflow mass rate is 
M'°"f=190 Mq yr“', consistent with the values observed for tar¬ 
gets at similar bolometric luminosities ( jCarniani et al.| 2015) 1. 
These equations assume a simplified model where the wind oc¬ 
curs in a conical region uniformly filled with outflowing ionised 
clouds. The values, not corrected for the extinction and regarding 
only the ionised component of the outflow (but see also the other 
conservative conditions in |Cano-Dfaz et al.|2012| and the discus¬ 
sion in |Perna et al.|2015] section 6.1), represent lower limits to 
the total outflow power. 

The kinetic power is 0.4% of the AGN bolometric lu¬ 
minosity also inferred from the SED fitting decomposition, in 
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Table 4: MIRO20581 outflow energetics 



basic asm 

min. asm 

max. asm 


(1) 

(2) 

(3) 

Pi'”' (erg s-‘) 

1.5x10^5 

3.8xl0« 

4.6x10'**^ 

M'Z (Mo yr-‘) 

2000 

48 

20000 

PriLboi 

0.04 

0.001 

0.4 

PZ'KLboilc) 

20 

0.5 

200 


Notes; Energetic values obtained from Eq.[ljand Eq. [fusing: basic as¬ 
sumptions (column 1), minimizing assumptions (column 2) and maxi¬ 
mizing assumptions (column 3) on n^, metallicity and extinction. 


iredictions of AGN feedback models 
The momentum flux, P™-Mvout 


rough agreeme nt with the 
(few % of Lboi, King|200^ 

is 2x10^® dyne, 2 times the radiative momentum flux from the 
central black hole, Lboilc. Its momentum boost, i.e. the ratio be¬ 
tween P‘™ and Lboilc is in agreement with those observed in X- 
ray winds (e.g., Tombesi et al.|20l3 1 rather than with the ratios 
associated with ionised and/or molecular outflows (e.g., |Cicone| 
|et al.|2014||Perna et al.|20l3| l. This discrepancy however, could 
be totally attributed to the fact that our estimate of the momen¬ 
tum flux represents a lower limit. In fact, correcting the [OIII] 
luminosity for an E(B-V)=1.3 calculated from the SED fitting 
decompositiot£] we obtain a mass outflow rate a factor of ~ 10 
larger (~ 2000 Mo/j/r) and, consequently, momentum flux and 
momentum boost in a more reasonable agreement with the re¬ 
sults reported in the literature (namely: a momentum boost 20 
times Lboijc, a kinetic power 4% Lboilc), and favouring an 
“energy-conserving” nature of the observed large-scale ionised 
outflow (see e.g. Carniani et al.|2015 and references therein). 
However, these results are based on few assumptions whose 
contribution can be remarkable. Estimates of the electron den¬ 
sity of outflowing region have been obtained in few low-z AGN 
and Ultra-Luminous Infrared Galaxies (e.g.,|Villar-Martm et al. 


2014 


et al. 


IRodriguez-Zaurm et al.|2013] l and high-z QSO (e.g., Perna 
|2015l Brusa et al. in prep), with values between l(F”and 


10^ cm Although the value we used in Eq. p] and Eq. i 


the one routinely used in literature (e.g., Liu et al. 2013[ Har¬ 


rison et al. 2014[ |Cresci et al. 2015| l, outflow energetics may 


be a factor of 10 lower if n^, =10^ cm is adopted. A further 
factor of 1/2 should be considered, taking into account the as¬ 
sumption on the metallicity, if metal-rich regions are present (see 
|Pema et^|2015| l. Einally, we considered a reddening value in 
between the E(B-V) estimated for the AGN and galaxy compo¬ 
nents (see Table |^. Considering the AGN (galaxy) reddening 
E(B-V) estimate instead of the average, the [OIII]5007 flux and 
therefore all the energetics should be corrected of a factor of 100 
(5). We note however that all basic assumptions previously out¬ 
lined are in agreement with those adopted in similar studies in 
literature and, in the case of the reddening value, reasonable. Eor 
completeness, we report in Tableall the energetic values with 
ranges obtained minimizing and maximizing the Eq.[T]and Eq.[^ 
using all the possible alternative assumptions above described. 


^ We note however that our kinetic power estimate is related to the 
only ionised component, while the prediction interest all the outflow 
components (i.e., molecular, atomic and ionised components). 

’ This value has been obtained as rough estimate assuming for the out¬ 
flowing material a reddening value in between the E(B-V) estimated for 
the AGN and galaxy components (see Table [^. In fact, possible dif¬ 
ferential obscuration between nuclear and off-nuclear regions may be 
present. 
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6. Discussion 

We analyzed the NIR SINEONI spectra of four candidate ob¬ 
scured QSOs, selected from the COSMOS survey, on the basis of 
red mid-infrared-to-optical and optical-to-near-infrared colours. 

Broad profiles in the [OIII] and Ha lines with EHWM > 550 
km/s (OC component) which are commonly used as signposts 
of outflows have been detected only in one source, MIRO20581. 
We can not exclude the presence of faint OC components in the 
other sources, given the low quality of the spectra and the lack 
of the [OIII]5007 emission, that is a better optical tracer for out¬ 
flows. Overall, the integrated spectra of the other three sources 
have low S/N (with emission lines detected at ~ 3 cr), and faint 
outflow components in the Ha-i-[NII] could still be present. Al¬ 
ternatively, the high obscuration of the sample, proofed the SED 
decomposition, the X-ray spectra, and the X-ray-to-mid-IR ra¬ 
tios (see Section [2T| l, might suggest that they are instead in the 
rapid black hole growth phase, when the scaling relations be¬ 
tween host galaxies and black hole properties are not yet estab¬ 
lished and winds have not been launched yet ( |King|2005] l. If it is 
the case, MIRO20581 is different from these other three sources. 
We investigate this scenario in the following section. 


6.1. Comparison with iiterature 


In this section we compare the properties of our MIRO targets 
with those of a sample of AGN-dominated DOGs reported in the 
literature, and discuss their main similarities/differences. Several 
studies have shown that DOGs with large f 2 Anm flux (> 1 mJy) 
exhibit higher AGN activity, higher concentration and smaller 
physical size. Vice versa, DOGs with lower f 2 Aum exhibit higher 


SE activity and larger physical size (Melbourne et al. 201 1[ 
Riguccini et al.|20T5 l. Although it is shown that the rest-frame 
optical morphologies of the most luminous DOGs have little sign 
of ongoing mergers ( [Melbourne et al.|2008 2009 1 Bussmann et 
al.|2009| l, there are also indications of non-regular gas kinemat¬ 
ics in their host galaxies ( [Melbourne et al.|201 Ij l and, moreover, 
several arguments suggest that they could be post-merger prod¬ 
ucts of gas-rich mergers (see e.g. Melbourne et al. 2009| l. Merger 
simulations were able to reproduce their colours and luminosi¬ 
ties and indicate an infrared-to-optical drop as gas consump¬ 
tion and AGN-driven wind terminate both SE and BH growth 
([Narayanan et al.|2010|l. 


Brand et al. ( 2007[ l studied NIR Keck spectra of a sample of 
ten AGN-dominated DOGs selected in the 9 deg^ NOAO Deep 
Wide-Eield Survey Bootes held Pannuzi & Dey[[l999[ ). Most 
of their spectra have low S/N and it is not possible to rule out 
the presence of outflows in the observed emission lines, but a 
prominent broad (EWHM ~ 1600 km/s) [OIII]5007 profile was 
detected in one target, SST24 J1428H-34. Ha or Wfi BLR have 
been found in 70% of them. [Melbourne et al.[ ( [201 l| l, instead 
presented a sample of four AGN-dominated DOGs with high 
spatial resolution Keck OSIRIS integral held spectroscopy. The 
sources were selected also in the Bootes held, and the main se¬ 
lection criteria was the strong Ha BLR detection in available 
NIR spectroscopic observations (e.g., from the same [Brand et al.[ 
2007 sample). They found that the BH masses of their sample 
are small for their host galaxy luminosities when compared with 
z ~ 2 and local unobscured AGNs (see their Section 5.1). Indeed, 
they did not And any evidence of outflows in the hosts, accord¬ 
ingly with the above-cited predictions of [King[ ( 2005[ l. Th ey also 


reported SERs <100 Mq yr ‘ for all targets. Einally, the 


Brusa 


et al. (20151 sample has been selected on the basis of red R-K 
colours (with a cut at K^b < 19), and high X/O ratio. The lat- 
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Fig. 8: (a): MIRO20581 contour plot of the median SINFONI data-cube over the entire HK wavelength range (green scale, in steps 
of Icr starting from 3cr), with over-imposed the Ks COSMOS (cyan; starting from 3cr, in steps of 3cr) and ACS HST contours 
(magenta; starting from 3cr, in steps of 3cr). (b): [OIII]5007 channel map obtained integrating the continuum-subtracted SINFONI 
data-cube on the total line profile (4973-5024A), of the region selected in panel (a). The contours levels are in steps of Icr (starting 
from 3cr) and show the [OIII] emission coming from the core (5000-5014A) and from the bluer (4973-5000A) part of the line 
profile, (c): [OIII] line profiles obtained from integrated spectra over the spaxels selected in the two regions ”A” and ”B” in the 
panel (b). Vertical dashed lines represent the systemic velocity obtained from the 1” integrated spectrum. 


ter selection criterion is roughly equivalent to the high MIR/O 
ratio (see [Fiore et al.|[7008| l. To confirm this we added the two 
brightest sources (XID2028 and XID5321; Perna et al.|2015j l in 
the [Brusa et al.j ( [2015 1 sample in Figure (with blue star and 
square), as representative of the entire sample. Evidence of out¬ 
flows have been found in 75% of object, and BLR emission have 
been found with similar percentage. 

Overall, all these sources have similar colours and MIR/O 
ratios to those of our MIRO targets (see Figure [T| cya n stars 
and squares represent the sources in [Melbourne et al.j p011| l 
and [Brand et al. ( 2007| l samples with K-band measurement). 
While BLR emission may be present only in 1 out 4 of our tar¬ 
gets (MIRO20581; see Secti on [5 1, a large number o f [Brusa et| 
[^ ( [2015) , [Melbourne et al.| ( |201 ]) and [Brand et al.| ( [2007| l ob-^ 
scured QSOs exhibit BLR emission. However, as already men¬ 
tioned, the Melbourne et al. p011| l targets were pre-selected to 
have a strong Ha detection and, in general, all the targets in 
these three comparative samples have strong K-band emission 
(Kab 20): the Brusa et al. (2015|l sample has been selected 


requiring a < 19; the Bootes field from which the Mel 
[bourne et al.| ( |201 1) and [Brand et al.[ ( [?007[ l DOGs were selected, 
has considerably shallower NIR observations than the COSMOS 
field (K/ib < 23, McCracken et al.|2010 < 20.8, Dey et al. 

[2008[ respectively). On the contrary, for the observations pro¬ 
posed in this work, we did not impose any flux threshold. Hence, 
it seems that the presence of BLR emission is related to the K- 
band flux. This offers a possible interpretation. All these targets 
are AGN-dominated DOGs at z~ 2; therefore, an higher K-band 
flux may correspond to an higher rest-frame optical AGN con¬ 
tinuum emission, that is proper of AGN showing BLR emission. 
Vice versa our sample, with Kab Z 20, is dominated by rest- 
frame optical host-galaxy continuum emission (see Ligure |^. 
The X-ray absorption is fully consistent with type 2 (1.9 in the 
case of MIRO20581) classification obtained by the SINLONI 
spectra. 


The only object in which we detected the outflow, 
MIRO201581, stands out with respect to the other targets. Al¬ 
though sharing the same X-ray luminosities of the [Brusa et al.[ 


( [2015| l targets, the X-ray spectrum of MIRO20581 shows high 
column density (N/y 7 x 10^^ cm“^), lar ger than those ob ¬ 


® cm-2. 


Perna et al. 


20151. 


served for the X/O targets (N/y 
This difference may be attributed mainly to line of sights effects 
which intersect a larger portion of the torus in MIRO20581 with 
respect to the Perna et al. ( [2015| l targets, and would also explain 
the fact that in MIRO20581 we do not detect a dominant BLR 
component in the Ha. Alternatively, the high extinction seen in 
the SED in both the AGN and host-galaxy components may be 
related to large scales obscurations. 


Assuming that the Ha-i-[NII] complex of MIRO20581 is 
composed of NC, OC and a BLR component, its black hole 
mass can be estimated using the Bongiorno et al. (20141 for¬ 
mula, assuming a LWHM of the Ha BLR emission of 3000 km/s 
(see Section[^, log(MB///Mo)=8.4. Considering the stellar mass 
estimated by SED fitting decomposition (log(M«/Mo)=l 1), we 
measure a central black hole to stellar mass ratio of the host of 
~0.002, comparable to the ratios of unobscured QSO at the same 
redshift ( Bongiorno et al.[[2014[ Merloni et al.|2010 l. Hence, it 
seems that outflows are present only in sources for which the 
Mbh - M, relation has been fixed, in agreement with the predic¬ 
tions of King (20051. 


Lrom the comparison between hboi (Table[^ and the Edding¬ 
ton luminosity associated to the M^h, we also infer an high Ed¬ 
dington ratio (LboilLEdd ~ !)• All these properties (high Edding¬ 
ton ratio, high extinction, M^h/M, ~ 0.002) point towards the 
interpretation that MIRO20581 could be associated with the be¬ 
ginning of the blow-out phase. On the other hand, its SLR and 
starbustiness are low (but still within +0.6 dex of the MS; see Ta¬ 
ble |g, indicating perhaps an advanced state of the ongoing pro¬ 
cess of negative feedback (see also [Balmaverde et al.|2015] l. We 
underline that its LIR emission is not well constrained (see Lig¬ 
ure [^; hence we calculated as upper limit a SLR = 132 Mo/yr, 
treating the LIR upper limits as real detections. An higher SLR, 
and consequently an higher sSLR, may be more consistent with 
the beginning of a blow-out phase, in which the effects of feed¬ 
back are still marginal and the SL is still occuri'ing. 
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In the [Narayanan et aT1 ( |2010) l simulations, AGN-dominated 
DOGs appear after the peak of the star formation, and therefore 
not necessarily associated with Starburst phase, given that the 
time scales of the two processes are not the same. Indeed, results 
recently presented in jRiguccini et al.] ( 2015| l showed that 50% 
of AGN-dominated DOGs detected by Herschel display sSFRs 
that place them in or above the MS, while the remaining 50% 
are below the MS, indicating perphaps an ongoing quenching 
of the star formation due to the AGN activity. In this scenario, 
small SFRs observed in dust-obscured massive main sequence 
galaxies hosting AGNs such those of MIRO10561, MIRO28704, 
MIR018744 and those presented in Melbourne et al. (2011 1 may 
be associated with systems which are still actively growing their 
black holes. 


7. Summary 

The main results from the multiwavelength analysis and the SIN- 
FONI data on a small sample of mid infraread bright, red quasars 
presented in this paper are summarised below; 


- All the sources but MIRO28704, selected from the 24;um 
Spitzer MIPS survey as candidate obscured QSO, are con¬ 
firmed highly (Nh ~ 2.5-7x 10^^ cm“^) obscured AGN from 
our detailed analysis. For MIRO28704 we found only indi¬ 
cations of a CT nature (See Section 2.1); 

- We successfully provided a spectroscopic redshift for three 
objects for which we had only a photometric estimate (see 
Section 4); 

- We revealed the presence of a powerful ionised outflow ex¬ 
tended out to ~ 4.8 kpc in only one source, MIRO20581 at z 
= 2.45. The large velocity (1600 km/s) and outflow mass rate 
(2000 Mo/j/r) for MIRO20581 are not sustainable by star 
formation. The energetics of the outflow are consistent with 
an energy-conserving mechanism (Section 5); the inferred 
Eddington ratio ~ 1), together with its highly obscured 
nature, point towards the interpretation that this source may 
be caught in the blow-out phase. 

- We collected several arguments that point towards the fact 

that luminous > 1 mJy) AGN-dominated DOGs may 

be objects in the transition phase between the post-merger 
starburst and unobscured QSO phases. The occurrence of 
outflows seems to be associated to the end of the rapid BH 
growth, when the - M, relation has been already estab¬ 
lished (See Section 6.1); 

- The efficiency of this mid-infrared-to-optical and optical-to- 
near-infrared colours selection criteria in detecting objects in 
the blow-out phase may be lower when compared to other se¬ 
lection criteria (e.g., [Brusa et al.|2015||Glikman et al.|2007) l. 
This may be due to the fact that higher obscuration could be 
associated to a still ongoing process of BH growth, preced¬ 
ing the blow-out phase (see Section 6). In order to confirm 
this statement, however, higher S/N spectra are required. Be¬ 
ing MIRO20581 the only target with outflow and an X-ray 
luminosity of 10“*^ erg s“’ (see Section 2), we suggest that 
sources in the blow-out phase can be most efficiently iso¬ 
lated from shallow X-ray surveys rather than solely on their 
high MIR/O colours (see alsojBrusa et al.|20r5]l. 


In the recent years, several detections of AGN-driven out¬ 
flows on kpc scale, probed by ionised gas kinematics, have been 
reported in the literature (see Carniani et al.|2015 for an updated 
compilation), and MIRO20581 is the last addition to this sparse 
and dishomogeneously assembled sample. Moreover, definitive 
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evidences of the the impact of the detected outflows on the host 
galaxies are still missing. In order to quantify how common are 
AGN driven outflows, and what is the impact of the wind in the 
hosts, NIR IFU AO-assisted observations of large (e.g. several 
tens) and homogeneously selected (e.g. from X-rays) samples of 
AGN are needed. We will address these issues in the near future 
through a SINFONI Large Program, “SUPER” (SINEONI Sur¬ 
vey for Unveiling the Physics and the Effect of Radiative feed- 
bacK, PI; V. Mainieri, ) which will target the first statistically 
sound sample (~ 40 AGN and QSOs drawn from the COSMOS, 
CDES and SDSS surveys) over four order of magnitudes in bolo- 
metric luminosities, and spanning all possible AGN (e.g. N//, 
Eddington ratio) and hosts (e.g. Starburstiness) properties. 
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